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The paramagnetic phase of the one-dimensional Kondo lattice model is investigated for electron 
densities below half-filling using a new mean-field approach. The physical parameters that govern 
this phase are identified to be the spin-flip processes of both the localized and itinerant spins. A 
new nonmagnetic-quantum state, where the local magnetization is a rotating vector with a nonzero 
average length, is proposed in order to describe this phase, this state does not break SU(2) sym- 
metry in agreement with Mermin- Wagner theorem. The line boundary between this phase and 
the ferromagnetic phase is calculated in the coupling-density phase diagram. Also, expressions are 
calculated for the velocities of the conduction electrons excitations, and heat capacity and entropy 
versus temperature are analyzed. Good agreement with many of the available numerical data is 
achieved. 
PACS numbers: 75.30.Mb 

I. INTRODUCTION 

Strong correlations between electronic degrees of free- 
dom in several materials, like high-temperature super- 
conductors, heavy-fermion materials, etc., are responsi- 
ble, or are at least thought to be responsible, for the 
occurrence of various quantum phenomena such as quan- 
tum magnetism, superconductivity, non-conventional 
metallic phases, etc. Heavy-fermion materials may be 
modeled using the canonical Kondo-lattice model for 
which intensive studies have been reported in the last 
years about its one-dimensional version. In this case, 
the zero-temperature phase diagram in terms of the 
conduction-electrons density n below half-filling and the 
Kondo exchange coupling constant Jk between the lo- 
calized and itinerant spins has been established mainly 
using numerical methods. Erlij For strong exchange cou- 
plings, the ground state is an unsaturated ferromagnet, 
whereas for weak couplings, it is a paramagnetic phase. 
At exactly half-filling, the phase is a spin-liquid insulator 
for any Kondo coupling, Fig. 1 (Ref. 1 contains detailed 
discussion of this phase diagram) . The ferromagnetic and 
the spin-liquid phases are now well understood because 
these are characterized by order parameters which are the 
magnetization for the ferromagnetic phase, and the_spin 
and charge energy gaps for the spin-liquid phase.El For 
the paramagnetic phase however, no physical parameter 
has been proposed so far in order to describe it. One of 
the purposes of this work is to identify such a parameter 
(or parameters). Moreover, the question concerning the 
crossover from the high-temperarture regime to the low- 
temperature region needs to be addressed. In the high- 
temperature regime, the conduction electrons and local- 
ized spins behave as being almost independent due to the 
fact that thermal fluctuations wash out any characteris- 
tic energy related to the Kondo exchange coupling. In 
the low-temperature regime, correlation effects become 
strong enough to dominate the physical properties of the 
Kondo lattice model. The paramagnetic phase was pro- 



posed to be a Luttinger liquid at zero temperature.liilul 
If we were to assume that this true, then a crossover 
would take place between the high-temperature metallic 
and the low temperature Luttinger-liquid phases. 

Our findings can be summarized as follows. The para- 
magnetic phase is found to be governed by the quantum 
dynamics (spin-flip processes) of both the localized and 
itinerant spins. We pave the way to a new quantum state 
in which the local magnetization is finite in the xy-plane 
if measured in a rotating reference frame, but vanishes 
once averaged over the angle of rotation, hence ensuring 
that Mermin- Wagner theorem is not violated; the mag- 
netization along the z-axis being equal to zero. This is 
an example of a rotating order prameter with a finite 
length but a phase angle assuming any value between 
and 2ir. We should mention that interpreting our find- 
ings in terms of the Luttinger liquid state is not a simple 
matter. We are still investigating this question. Note 
however that quantitative comparison between our re- 
sults and the Luttinger-liquid description results is found 
to be fairly good (see below) . What remains to do is un- 
derstanding: why should the Fermi liquid picture give 
rise to the Luttinger liquid one at low temperature? 

In section II, a full description of the present approach 
is provided. The physical foundations upon which our 
theoretical calculations are based, are explained. We 
identify the physical parameters as being the averages 
of spin flip operators of the localized and itinerant elec- 
trons. Section III is devoted to showing some results and 
their discussions. The temperature dependence of the pa- 
rameters, crossover temperature, phase-diagram bound- 
ary between the ferromagnetic and paramagnetic phases, 
as well as some of thermodynamic functions are reported. 
Comparison with numerical data turns out to be very sat- 
isfactory. In section IV, conclusions are drawn about the 
validity of the present approach. 
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II. DESCRIPTION OF THE APPROACH 

A. Physical foundations 

In this work, the numerical data available in the litera- 
ture are considered as numerical experiments upon which 
the present approach is founded. 
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FIG. 1. The phase diagram obtained using Eq. 14 is shown 
and compared to numerical results. FM, ferromagnetic phase. 
PM, paramagnetic phase. 

The Kondo lattice Hamiltonian in one dimension is 
given by: 

h = -t Y,( c U c i+ho + H - c ) + j k s * • (*) 

Z,f7 i 

where and are respectively the localized and itin- 
erant spin operators. c\ a and Ci. a are the creation and 
annihilation operators at site i of an electron with the 
z-component of the spin being a = ±1/2. t is the 
hoping energy of conduction electrons. Because neither 
the Kondo singlets formation nor magnetic order oc- 
curs in the paramagnetic phase for densities below half- 
filling, a result that is strongly pointed out by numerical 
calculations,^ it is justified to choose the following canal 
for decoupling the interacting term of ([l]) where only the 
averages (s~) and (S~) and their complex conjugates are 
taken into account. These parameters represent spin- flip 
processes, and are therefore a good probe of the quantum 
dynamics of the system. We suppose that these processes 
dominate over all other processes which are the Kondo 
screening and magnetic ordering. We will see that our re- 
sults turn out to be consistent with this hypothesis. It is 
worth noting here that in order to avoid breaking SU(2) 
symmetry, the phase angles of these parameters will not 
be calculated in mean-field theory. We, rather, perform a 
summation over all possible values of these angles. This 
procedure ensures that the physical phase we obtain is 
nonmagnetic because it does not break SU(2) symmetry. 

B. Effective Hamiltonian 

Using this new canal, we get the following approximate 
expression for the interacting term of (pi): 



S, • 8i « -[ST A* + s rQ* + H.c. - 2Re(QiA*)], (2) 

where Q t = (Sf) = \Qi\e** and A, = («r) = \A t \e^ . 
Then, the moduli \Q%\ and \Ai\ are calculated in mean- 
field approximation, but summation over the phase an- 
gles ipi and fa a re performed to guarantee that the con- 
tinuous SU(2) symmetry remains unbroken. The to- 
tal Hamiltonian is averaged over fa and ipi, and the 
summation over the phases fa and ipi of the lowering 
spin operators lead to (S^)<p i = dfa \ Qi \ e l< ^ i = 

and (Sj - )^ = J 27r dip^Aile 1 ^ = where \Q t \ and fa 
on one hand, and \Ai\ and ipi on the other hand are 
considered to be independent variables. The minimiza- 
tion of the average of the magnetic energy Jx(S,Si) = 
</if cos(fa — ipi) imposes the constraint fa — ipi = ^ 

on the angles. For this reason, for example, the sum over 
these phases in the last term of (Q) , which involves A* Qi , 
is / difa-iPt-njlAtQ^e^-^Sifa-iPt-Tr) = -\AiQi\. 
The delta function implements the constraint fa— ipi = n. 

Using the second quantization form for the itinerant 
spins, the Hamiltonian ([!]) leads to the following effective 
Hamiltonian where averaging over the phase angles ipi 
and fa is done: 

2 i Jo 27T 

-t 2_j 4,<r c i+hcr + H.C 

~Y1 j - <l ) i) X 

i 

^{lAQile^-^Sifa -iPi-Tr) + C.c.}. (3) 

C. Unifrom rotating configuration 

In the rest of this work, we focus our attention on the 
uniform configuration, which is obtained for Qi = \Q\e lc ^ 
and Ai = \A\e 1 ^ . Q, A, <p and ip are considered to be 
site independent but angles vary between and 27r while 
satisfying the constraint <p — ip = it. This state realized 
in this way is called the uniform rotating configuration 
(URC). In this state, the vector parameters Qi = (S^~) 
and Ai — (s~) point in opposite directions while rotating 
at the same rate. The local and total magnetizations are 
equal to zero. 

In addition, we treat the up and down fermions as 
being different, and use the following canonical transfor- 
mation: 

c kl = e^ 2 (p k +a k )/V2 
c kl =e^ 2 (p k -<j k )/V2 (4) 
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in order to diagonalize the effective Hamiltonian (||) in 



the case of the URC. Here, = £\ c\ v e~ ir - k / VN with 



being the number of lattice sites. Next, we perform 
the following transformation to absorb the phase terms 
in the transformations (0): 



e^ /2 c fc , T , 



Cfc.l 



40/2 



(5) 



which is equivalent to making a rotation by angle cp about 
the z-axis for the x- and y-components of the itinerant 
spin operator, with the matrix of rotation given by 



cos< 
— sine 



sin<j 
cos< 



Due to the constraint <f> — ip = tt, rotating the itinerant- 
spins x and y components causes the rotation of the x and 
y components of the localized spins by angle ip — (p — tt 
about th z axis. The matrix of rotation in this case is 
given by: 



cos ip sin ip 
sin ip cos ip 



Because we sum over the angles <p and ip between and 
2-7T, the spin components are continuously rotating, and 
the Hamiltonian may be written in the rotating reference 
frame for which the x axis and y axis coincide with the 
rotating x and y components of the itinerant spins. 



D. Mean field Hamiltonian in the rotating reference 
frame 

The fact that the effective Hamiltonian is invariant un- 
der the above rotations by angles <p an d ip is consistent 
with the absence of SU(2) symmetry breaking. Here we 
perform such rotations before diagonalizing the simpli- 
fied Hamiltonian obtained in the URC approximation. 
This is equivalent to using the transformations 



e^ /2 c, T , 



(6) 



on both the itinerant and localized spins. The result is a 
much simpler expression for the Hamiltonian: 

k 

+aj k Y / s? + JkY, a Q ( ? ) 

i i 

where A and Q stand now for the magnitudes \A\ and 
\Q\ respectively. Here, E Pt<r = e(k) — \i ± \Q\Jk/2; 
e(k) = — 2icosfc is the tight-binding spectrum, and fi 
is the chemical potential. Note that we should keep in 
mind that (|7]) is obtained in a rotating frame as <j> and ip 
take values in the interval [0, 2tt{ subject to the constraint 



<p — ip = 7T. Therefore, it is inappropriate to interpret this 
Hamiltonian as that of conduction electrons and localized 
spins coupled to magnetic fields along the x-direction in 
a reference frame at rest. Now, we understand that the 
up and down spins could be treated differently as a con- 
sequence of the effective rotating magnetic fields AJk 
for the localised spins, and QJk for the itinerant spins. 
Note also that the energy spectrum of the itinerant elec- 
trons splits into two bands under the effect of the Kondo 
exchange interaction. 



E. Phase fluctuations 

Allowing for the summation over the phase angles 
means that we do not minimize free energy with respect 
to these angles. However, to satisfy the requirement 
that the magnetic energy is minimized we constrained 
their difference to be w. This leads us to questioning 
whether this minimum is stable againt fluctuations about 
the value <j> — ip = ir. To study the effect of these fluctu- 
ations, the following treatment is done. We replace the 
delta function in the last term of TL in Eq. (||) by the 
broader Gaussian distribution 
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-(0-^-7T) 2 /f 2 



to allow for other angles difference to contribute. If e — > 
the Gaussian distribution reduces to the Dirac distribu- 
tion. The last term in (0) takes the form: 



d(tP-iP) 



Jk 



-(0-V-7T) 2 /6 2 



C.C.}. 



Integration over e from — oo to +oo is undertaken, once 
the integration over <p — ip is done, to guarantee that all 
possible fluctuations of the phase angles are embodied in 
the present approach. The result in the URC is 



-aJ K J2\AQ\ 



which differs from the result — Jk Yli \AQ\, obtained 
without fluctuations, by the factor a. a s a number 
larger than 1 but smaller than which is the value 
obtained when the integrations are carried out to infinity 
in the integral on (p — ip. This reduces the values of the 
mean-field parameters (as one would expect) by a factor 
a ~ 1 without destroying the mean-field picture. Thus 
this means that the mean-field solution <p — ip = n is 
stable against phase fluctuations. 



III. RESULTS 
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A. Parameters of the KLM model 



B. Crossover temperature 



In the rest of this paper, we set a = 1, and seek some 
quantitative understanding of the present approach. To 
calculate the parameters Q and A, we use the self- 
consistent equations obtained by minimizating the free 
energy: 



F = - 



h E >»<i 



N/3 



} 



--ln{2cosh[/L4J K /2]} + J K QA. 

P 

with respect to A and Q. This yields: 
Q = -t&nh(pj K A/2) 



(8) 



(9) 



where the summation runs over the Brillouin zone of 
the conduction electrons, and f(x) = 1/(1 + e^ 35- ^) is 
the Fermi distribution factor. (3 — 1/ksT with ks is the 
Boltzmann constant, fi w e F with e F being the Fermi 
energy of the conduction electrons. 

As the paramagnetic phase occurs for Jk < 4t, we ex- 
pand f(E p ) — f{E a ) to first order in QJ K /2. We obtain: 



.4 



1 



df(e(k)) 
de(k) 



1 



K ! 



(10) 



which leads to 



The parameters A and Q decrease as temperature in- 
creases, and vanish at a temperature T c f. Indeed, ex- 
panding the tangent hyperbolic to third order in the 
vicinity of this critical temperature in Eq. ( pi] ) and using 
Eq. (y), we find: 



1/2 



A ~ A Q {1 - T/T cf }V\ Q ~ Q {1 - T/T cf } 

with Aq = 2v / 3fcsT' c // Jx, and Qq = V3/2, and the tem- 
perature 



Tcf = xJk/Mb- 



(12) 



Below T c f, the Kondo lattice system deviates from the 
system of independent localized spins and conduction 
electrons. As might be expected T c t is proportional to 
xJki a result that is reminiscent of the RKKY interac- 
tion which dominates over Kondo screening, although not 
leading to magnetic order. Above T c f, the magnitude of 
the magnitization is zero. Below, T c f, the magnitude be- 
comes finite, but averaging about its phase angles gives 
zero magnitization. It is for this reason that this change 
in behavior is not a true phase transition. The present 
appoach leads to a sharp change of regime at T c f. We 
do not exclude that if the change in the regime is rather 
a smooth crossover, the crossover temperature will be 
given by T x ~ T c /. T c f is called the crossover tempera- 
ture. Shibata and Tsuntsuguo reported the existence of 
a crossover temperature, but said that it is determined 
by the velocity of the excitations of the itinerant spins at 
zero temperarture. 



A= - X J K t&nh(f3AJ K /2) 



(11) 



when 



is used. Here x = — J2k 



is the uniform 



df{e(k)) 
■>k de(k) 

susceptibility of the free conduction electrons in one di- 
mension. 

At zero temperature, the non-zero solution is given by 
A = xJk/4 and Q = 1/2; X ~ 1>{e F ) = 1/2ttv f is the 
density of states of the free conduction electrons; vp = 
2isinfci? being their Fermi velocity. Slightly above zero 
temperature, we get for fc^T <C V{tp)j\j ^ 



A~ - x Jk< 
„ 1 



1 - 2e 



where again \ ~ ^^f) is evaluated at T = because 
x(T) does not deviate much from its zero-temperature 
value for temperature well below the Fermi temperature 
of the free conduction electrons. In the rest of this paper 
we will consider this value for X - 



C. Ground-state energy and the phase diagram 

The idea behind calculating the ground-state energy 
is to find how the correction (due to Jk) to the ground- 
state energy of the conduction electrons behaves. This 
will indicate the type of correlations that dominant. This 
correction is found to be given by the equation: 



AE, 



GS 



(\AE GS \/k B T cf w 1) 



(13) 



for Jk <C 4i. Interestingly, this correction is of order 
J K 1 1, and is again reminiscent of the RKKY interaction 
between localized spins. As Jk /t increases (but keeping 
Jif/i < 4), we found that the correction to the ground- 
state energy deviates from a J K /t law. Estimating the 
points where the change takes place for different values of 
Jk, we determined the line boundary, JK,c(n)/t versus 
the conduction electrons density n, separating the para- 
magnetic phase from the ferromagnetic phase at strong 
coupling JK/t- This line is found to be very well fitted 
by: 



JK,c(ji)/t ~ 2im sin(7rn) . 



(14) 
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This is consistent with the fact that the weak coupling 
regime is equivalent to Jx/t < ep/t. Linearizing the 
energy spectrum of the conduction electrons around kp, 
one finds the linear spectrum e(k) — ep — vp{k — kp), 
where vp is the Fermi velocity. For a particle with mo- 
mentum k {Ti — 1) and velocity vp, the energy is given 
by vpkp. Thus, within the linear approximation, it is 
natural to consider an effective Fermi energy given by 
ep = vpkp = Itkp sinfcp. This yields Jk < vpkp = 
2tim sin(7rn), which leads to (|f4|). 

Eq. (14) is in very good agreement with the avail- 
able exact numerical data, Fig. 1. For n — 0.35, Eq. 
(14|) yields Jk.c ~ 1-96, which is very close to 



Av = v n 



Up ~ ^r- cot(rar). 



(15) 



density- matrix-renormalization-group exact result, 2li 
The agreement is however less accurate for n = 1/6 whpe 
we obtained Jk.c ~ 0.5 compared to the result of Ref.t3, 
namely Jk.c ~ 1, obtained using quantum Monte Carlo 
simulation. The reason for this discordance may how- 
ever be attributed to the fact that Monte Carlo simula- 
tions are difficult to implement at very low temperatures, 
and to the fact that the present theory is mean-field like. 
Using the bosonization techniqueEj Honner and Gulacsi, 
reported J/f iC (n) « 2.5 sin(7m). But as we notice on fig- 
ure 1, this does not capture the general trend, while Eq. 
( |l4| ) does it fairly well. Honner and Gulacsi noted that a 
dependence as that I report here is also possible. 




FIG. 2. The difference between the spin and charge veloci- 
ties of the conduction electrons is displayed in terms of Jk/1- 



D. Excitation velocities 

For densities below the half-filled band, and in the pres- 
ence of the Kondo coupling, we have found that the con- 
duction electrons as charge and spin entities still form a 
gapless phase because the term ±\Q\Jk/2 in E P:(T does 
not open a gap in the energy spectrum. The elementary 
excitations are described by the operators pu and Ok- 
Their excitation velocities are given by v P:Cr = 2tsuik p<7 
where k Pia — nn =p J/f/8i sin(n7r) for small Jk- k p 
and k a are determined by the conditions E p = ep and 
E a = ep with ep = — 2t cos(kp). Thus, the p- and 
CT-velocities are given by v p ~ vf — Jk cot(mr)/4 and 
v a ~ vp + Jk cot(n7r)/4 for Jk < JK,c{n). This leads to 



It is not clear how we could relate these velocities to 
the spin and charge velocities of the conduction electrons 
obtained within the bosonization approach. In Fig. 2, we 
draw Av in terms, of Jk for n — 1/3, and compare it with 
Shibata et aZ.'alj density- matrix-renormalization-group 
results. Note that this comparison should be considered 
with caution as our results for v p and v a disagree with 
theirs. Only the difference compares well. We should also 
emphasize here that we do not pretend by any means 
that the present mean-field theory is equivalent to the 
bosonization technique which is used to obtain the spin 
and charge velocities when spin-charge separation takes 
place in a system of interacting fermions. 



E. Thermodynamic functions 

Finally, we would like to briefly analyze the entropy 
S = -dF/dT and heat capacity C = TdS/dT. In Fig. 
3, we display C, for n — 0.35 and Jn/t = 1.6 where 
T c f = 0.057i/fcs agrees very well with the crossover tem- 
perature predictfijd by Shibata and Tsunetsugu in their 
numerical worked. The entropy S is also reported on 
the same figure. It behaves as predicted in RefEI For 
T/T cf — > 0, the entropy S ~ 1.24T yields a slope in vecyi 
good agreement with the Tomonaga-Luttinger resultEll 



S = nT(v 



)/3 ~ 1.20T. A sharp maximum is 



found in C at T c f as a consequnece of the change in the 
regime from hight temperatute to temperatures below 
T c f. The linear behavior of entropy in terms of temper- 
ature is consistent with the gapless quasi-particle exci- 
tations, and is an indication of the metallic character of 
the paramagnetic state. 

2 
1.5 
1 

0.5 


0.1 0.2 0.3 0.4 

k B T/t 

FIG. 3. The specific heat and entropy are drawn vs T/t for 
Jk = 1.6t and n = 0.35. Comparison is made to the results 
of Ref. 20. 



IV. CONCLUSIONS 

In summary, the paramagnetic phase of the one- 
dimensional Kondo lattice model is investigated using 
a new mean-field approach. This approach is valid 
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away from half-filling where Kondo screening is negli- 
gible. We predict that this phase is described by a new 
nonmagnetic-quantum state characterized by two rotat- 
ing order parameters with nonzero magnitudes below a 
crossover temperature. We calculated the ground-state 
line boundary in the density-coupling phase diagram, the 
elementary-excitations velocities of the conduction elec- 
trons, heat capacity, and entropy. The crossover tem- 
perature separating the low-temperature paramagnetic 
quantum phase from the normal high-temperature phase 
is evaluated as a function of the Kondo coupling con- 
stant and the density of states of the conduction elec- 
trons. Overall, our results agree very well with many 
of the available numerical data well below half-filling, 
and are compatible with the Luttinger-liquid picture as 
put forward by several authors, (see Ref.u and references 
therein) . Note that it is not yet clear how we could inter- 
pret the results of our approach in terms of the Luttinger- 
liquid description. For densities close to half-filling, resid- 
ual effects due to Kondo screening have to be taken into 
account. It is then natural that, quantitatively, our ap- 
proach is less accurate in this limit. This is what we 
come tctface when we compare Shibata and Tsunetsugu's 
resultsE3 with ours for the entropy S for example. The 
results start to agree well only for conduction-electrons 
densities below 0.35. Finally, the RKKY oscillations con- 
stitute an issue that needs to be addressed within the 
present approach. One possible avenue is to allow for 
the parameters A Q to be fc-dependent. However, even 
with these two limitations, the present approach is very 
promising because it is very simple to use, and the calcu- 
lations can, to a very large extent, be done analytically. 
And most importantly, it leads to many very satisfactory 
physical results. 

The author would like to thank Prof. P. Fulde for his 
comments on the manuscript. 
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